Implications of the AMS-02 positron fraction in cosmic rays 
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The AMS-02 collaboration has just released its first result of the cosmic positron fraction e+/(e^ + 
e"*") with high precision up to ~ 350 GeV. The AMS-02 result shows the same trend with the previous 
PAMELA result, which requires extra electron/positron sources on top of the conventional cosmic 
ray background, either from astrophysical sources or from dark matter annihilation/decay. In this 
paper we try to figure out the nature of the extra sources by fitting to the AMS-02 e"'"/(e~ -|- e"*") 
data, as well as the electron and proton spectra by PAMELA and the (e~ -I- e^) spectrum by Fermi 
and HESS. We adopt the GALPROP package to calculate the propagation of the Galactic cosmic 
rays and the Markov Chain Monte Carlo sampler to do the fit. We find that the AMS-02 data 
have implied essential difference from the PAMELA data. There is tension between the AMS-02 
e"'"/(e~ -I- e"*") data and the Fermi/HESS (e~ -I- e"*") spectrum, that the AMS-02 data requires less 
contribution from the extra sources than Fermi/HESS. Then we redo the fit without including the 
Fermi/HESS data. In this case both the pulsars and dark matter annihilation/decay can explain 
the AMS-02 data. The pulsar scenario has a soft inject spectrum with the power-law index at ~ 2, 
while the dark matter scenario needs t'^t~ final state with mass ~ 600 GeV and a boost factor 
~200. 

PACS numbers: 95.35.-|-d,96.50.S- 



I. INTRODUCTION 

The Alpha Magnetic Spectrometer (AMS-02) was 
launched in May 2011. After nearly two years oper- 
ation and analysis, the AMS-02 collaboration has re- 
leased its first physical result, i.e. the positron fraction 
e+/(e^ + e+) in cosmic rays (CRs) [l[. The data show 
very high precision with the energy range from ~ 0.5 
GeV to - 350 GeV. The fraction rises above - 8 GeV 
up to the energy end at ^ 350 GeV, which is consistent 
with the previous PAMELA result of the cosmic positron 
fraction [2J,0] . The result is NOT consistent with the con- 
ventional CR expectation [j] . A great number of works 
have dedicated to explaining the PAMELA result, either 
by astrophysical sources J5l-[l0| or by dark matter (DM) 

[nHH. 

To determine the parameters of the CR background 
and the nature of the extra sources we have to consider all 
the relevant results. The available results at present in- 
clude the pure electron spectrum measured by PAMELA 
[3 and the antiproton flux and ratio p/p [13, [l3l by 
PAMELA. Another important result from PAMELA is 
the precise measurement of the proton spectrum [T^, 
which determines the secondary positron spectrum by 
collision with the interstellar medium (ISM) when prop- 
agating in the Galaxy. There are also precise mea- 
surements of the total electron and positron (e" -I- e^) 
spectrum, by the Fermi-LAT collaboration |20l . [21| and 
ATIC collaboration [22] ■ The ground- based atmospheric 
Cerenkov telescopes HESS also gives measurement of the 
total e^ spectrum up to higher energies |23l.[24|. 

Fitting to the PAMELA positron fraction data [g, 
PAMELA electron spectrum and Fermi/HESS total e^ 



spectrum shows that both the astrophysical source, such 
as pulsars, and the DM scenarios can give a good expla- 
nation to the data [J|. It is hard to discriminate the two 
scenarios with the CR spectra mentioned above. As the 
AMS-02 data show much higher precision and wider en- 
ergy extension, especially it shows softer behavior than 
the PAMELA 2008 result, it is necessary to re-examine 
the previous conclusion with all of the newly available 
data. In this work we have done such a global fitting 
to all the relevant data, including AMS-02 e"'"/(e~ -t-e"*"), 
Fermi and HESS (e" -|-e+) total spectrum and PAMELA 
proton and electron spectra. 

We adopt the CosRayMC code, which embeds the 
CR propagation code in the Markov Chain Monte Carlo 
(MCMC) sampler and enables efficient survey of the 
high-dimensional parameter space [J, [23I ■ The CR prop- 
agation is treated by the GALPROP package [2a|. The 
CR transportation process in the Galaxy is character- 
ized by the secondary particles. Therefore the secondary- 
to-primary ratio, such as B/C, (Sc-|-Ti-|-V)/Fe, and the 
unstable-to-stable ratio of secondary particles, such as 
^°Be/^Be, ^^Al/^^Al are often used to determine the 
propagation parameters [26l - l29l | . In this work we fix the 
CR propagation parameters to the values which give the 
best fitting to the currently available B/C and ^°Be/^Be 
data with the MCMC method. The fitting process will 
be reported separately in [3(j |. 

We then fit the parameters of the electrons and 
positrons to the relevant data, both from the CR back- 
ground and the extra sources. In the work we have 
considered the continuously distributed pulsars and DM 
annihilation/decay as two typical scenarios of the extra 
positron/electron sources. The injection parameters of 



primary electrons are free parameters to be fitted. The 
cosmic positrons include secondaries from CR interac- 
tion with the ISM and the primary ones from the extra 
sources. The secondary positron spectrum is determined 
by the spectrum of cosmic protons (including Helium and 
a few heavier nuclei) and the primary positron spectrum 
is determined by the nature of the extra sources. In prin- 
ciple the injection parameters of the protons can be fitted 
independently and then be employed to calculate the sec- 
ondary positrons. However, since both the proton and 
electron spectra are measured by PAMELA at almost 
the same time, they are modulated by the solar activ- 
ity with a similar magnitude. Therefore, in this work we 
have adopted two ways to fit the cosmic proton injection 
parameters: either by fitting the proton spectrum inde- 
pendently or by fitting the proton and electron spectra 
simultaneously. In the first case we enable the protons 
and electrons/positrons to have different solar modula- 
tion parameters, while in the second case we assume they 
share the same modulation parameter. 

By the global fitting we find it seems not easy to fit all 
the relevant data simultaneously. The AMS-02 data and 
the Fermi/HESS data might have a tension between each 
other. The AMS-02 requires less contribution from the 
extra source while Fermi data requires more. We then do 
the fit without including the Fermi/HESS data. We find 
in such case both the pulsar and DM scenarios can give 
good description to the data. 

This paper is organized as follows. We give a brief 
introduce of the propagation of Galactic CRs in Sec. [Hi 
The experimental data and fitting method are described 
in Sec. IIIII The results are presented in Sec. IIVI In Sec. 
|V]we give discussion about the fitting results, and finally 
a brief summary is given in Sec. I VII 



II. COSMIC RAY PROPAGATION 



with the spatial diffusion coefficient D^x as [32 



^PP^--- SS{A-S^){A-S)w' 



(2) 



where va is the Alfven speed, w is the ratio of magneto- 
hydrodynamic wave energy density to the magnetic field 
energy density, which characterizes the level of turbu- 
lence. In the usual way we take w to be 1 and use the 
Alfven speed va to describe the reacceleration (s^. The 
CRs propagate in an extended halo with characteristic 
height Zfi, beyond which free escape of CRs is assumed. 
Thus the major propagation parameters include Dq, S, 
VA, Vc and Zh- 

There are publicly available numerical codes to com- 
pute the CRs propagation in the Galaxy, such as GAL- 
PROPi [111 and DRAGON^ ^. We have embedded the 
GALPROP package with the MCMC sampler and fitted 
to the data to constrain the propagation parameters |30[ . 

Recently several groups employed the MCMC tech- 
nique to fit the CR propagation parameters J34| - |36| . Us- 
ing the currently available data of B/C and ^°Be/^Be, we 
did an independent MCMC fit to the propagation param- 
eters [301 • We find the reacceleration model gives quite 
good description to the present data, while the convec- 
tion model gives worse fitting. Therefore we adopt the 
reacceleration (i.e., Vc = 0) scenario as the starting point 
of the present study. The values of the propagation pa- 
rameters we adopted in the work are listed in Table |T[ as 
given in [30|. 

TABLE 1: Propagation parameters taken in the work 
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The propagation of charged CRs in the Galaxy is de- 
scribed by the diffusive equation [31| 
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where ip is the density of CR particles per unit momen- 
tum interval, Q(x, p) is the source term, D^x is the spatial 
diffusion coefficient, Vc is the convection velocity, Dpp 
is the diffusion coefficient in momentum space used to 
describe the reacceleration process, p = dp/dt is the mo- 
mentum loss rate, tj and r^ are time scales for fragmen- 
tation and radioactive decay respectively. Dxx is usually 
assumed to be only rigidity dependent and has a power- 
law form Dxx = DqJ3{R/ RqY , with 5 reflecting the prop- 
erty of the interstellar medium (ISM) turbulence. The 
reacceleration is described by the diffusion in momentum 
space. The momentum diffusion coefficient Dpp relates 



III. FITTING DESCRIPTION 

A. Model 

In the conventional CR model the positrons are 
secondary particles generated by the interaction be- 
tween the primary nuclei and the ISM. Wide studies of 
the PAMELA result show that, besides the secondary 
positrons, extra sources of positron/electron are nec- 
essary to explain the positron fraction e+/(e^ + e+) 
3. The present model to fit the AMS-02 data also in- 
cludes the traditional CR background (primary nuclei 



^ http://galprop.stanford.edu/ 

^ http://www.desy.de/maccionc/DRAGON/ 



and electrons, and the secondary positrons) and the extra 
positron/electron sources. 

The injection spectra of the primary protons (heavier 
nuclei are less relevant) and electrons are assumed to be 
broken power-law functions with respect to the momen- 
tum p 



lip) °^ i^^] 
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where v\ and vi are the spectral indices below and above 
the break momentum pbr- In the following we use v\ and 
vi to represent the proton spectrum and use 71 and 72 
to represent the primary electron spectrum. The propa- 
gated fluxes of protons and electrons are then normalized 
to factors A^, and v4e to get the absolute fluxes. 

The spatial distribution of the primary CR particles is 
adopted to be the supernova remnants (SNR) like distri- 
bution 
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where Rq =8.5 kpc is the distance of solar system from 
the Galactic center, Zg pa 0.2 kpc is the characteristic 
height of the Galactic disk. We employ a = 1.25 and 
b — 3.56 following [Sa]. 

We employ the parameterization given in [37| to calcu- 
late the secondary electron/positron spectra from pp col- 
lision. Another free factor Ce+ to normalize the secondary 
positron/electron flux is included in the fitting, which 
represents the possible uncertainties from the hadronic 
interactions, propagation models, the ISM density distri- 
butions, and the nuclear enhancement factor from heavy 
elements. 

About the extra electron/positron sources, we study 
two popular kinds of sources: the astrophysical sources 
such as pulsars and the DM annihilation scenario. The 
injection spectrum of pulsars is parameterized as a power- 
law function with an exponential cutoff 



q{p) = ApsrP "expi-p/pc), 



(5) 



where Apsr is the normalization factor, a is the spectral 
index and Pc is the cutoff momentum. The pulsar popu- 
lation is taken to be a continuous form with distribution 
function ((4]), but with different parameters a = 2.35 and 
b = 5.56, given in [3g. The effect of nearby isolate pul- 
sars (e.g., [ll-Q, see also the discussion in Sec. V) is not 
covered in the present study. 

For the DM annihilation scenario, the density profile of 
the Milky Way halo is adopted to be the Navarro-Frenk- 
White (NFW) distribution H 
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with Vg = 20 kpc and Ps — 0.26 GcV cm"'^. Such a 
value of Ps will correspond to a local density of 0.3 GcV 
cm""^. For higher values of the local density as revealed 



by several recent studies |40l - l42| . the annihilation cross 
section will be different by a constant factor. Since the 
measurement of CR antiprotons by PAMELA [13, [Hi 
constrain the hadronic annihilation channels strongly 
prl Il2l . |43| . |44| , we will focus on the leptonic annihila- 
tion channels here. The positron spectrum from the DM 
annihilation products is calculated using the PYTHIA 
simulation package [45 1. 



The CRs at low energy (typically with rigidity below 
~ 30 GV) are affected by the solar environment when 
entering the solar system, known as solar modulation ef- 
fect. The force field approximation is often employed 
to describe the solar modulation effect [4^, which has 
only one single parameter — the modulation potential (j). 
However, the low energy data about the positron fraction 
measured by PAMELA and AMS-02 may imply that the 
simple force field approximation is not enough to explain 
the data, and the charge- sign dependent modulation ef- 
fect is necessary |47H49| . Therefore, in order to avoid 
possible inconsistency of the low energy behavior we do 
not include the AMS-02 data below 5 GeV in our fit. The 
solar modulation effect above 5 GeV is negligible [4^ [50] . 

In summary the full parameter space investigated in 
this work is 



P = 



' {Ap,Vi,V2,pl^}, 
Me, 71, 72, Kr}' 

{A.psr,a,pc} or {m 

.{Ce+,0}, 



X' ' 



^)}, 



bkg protons, 
bkg electrons, 
exotic sources, 
others. 



(7) 



B. Data 



In this study the data to be fitted include the latest 
positron fraction by AMS-02 [l[, the electron spectrum by 
PAMELA (ig . the total electron and positron spectra by 
Fermi [23,[2l| and HESS [2i,[2i|. As for the proton spec- 
trum which will affect the secondary positron production, 
we choose two ways to deal with it as described in Sec. L 
to fix it to the PAMELA measurement [1^ or to include 
the injection parameters in the fit. 

Note that the CR spectral hardening at ^ 200 GV 
reported by ATIC [Hj, CREAM H^ and PAMELA [H 
implies that single power-law can not fully describe the 
high energy spectra of the CR nuclei above ~ 10 GV. 
Therefore we only take the PAMELA proton data below 
150 GeV in the fit. To better fit the high energy part 
(and the CREAM data), we need a further break or a 
curved injection spectrum of the protons [53[. 

Since the Fermi and HESS data may have larger sys- 
tematic uncertainties we also investigate the case without 
the Fermi/HESS data. The definition of the fittings are 
compiled in Table HIl 



TABLE II: Definition of fitting 



I-a 
Il-a 



AMS e+/e± + PAMELA e 
AMS e+/e± + PAMELA e 



Fermi/HESS e^ 



AMS e+/e± + PAMELA e^ + Fcrmi/HESS e± + PAMELA p 
AMS e+/e^ + PAMELA e~ + PAMELA p 



I-b 
Il-b 



IV. RESULTS 

A. Fixing the proton spectrum 

We first fit the proton spectrum independently with 
the PAMELA (and CREAM) data. Here we add the high 
energy CREAM data [52| in the fit to give a description 
of the proton behavior in a wider energy range. The best 
fitting injection parameters of protons are: vi = 1.79, 
1/2 = 2.36, p^j. = 11.7 GeV and the solar modulation 
potential (f> = 470 MV. The propagated proton spectrum 
for the best fitting parameters is shown in Fig. [T] 
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FIG. 1: Proton spectrum derived through fitting PAMELA 
data. References of the proton data: AMS [111 , BESS [111 , 
ATIC2 [H, PAMELA ^ and CREAM ^. 

We then run the fits I-a, Il-a to derive the parameters 
of CR electrons/positrons with the best fitting proton 
spectrum. The resulting positron fraction and electron 
spectrum for the best fitting parameters of each fit, for 
both the pulsar and DM scenarios, arc shown in Figs. [2] 
- m For each figure, the panels from left to right cor- 
respond to the fits I-a and Il-a respectively. The best 
fitting parameters, mean values and their la uncertain- 
ties are compiled in Tables IIIII - |Vl and the best fitting 
X^ over the number of degree of freedom (dof ) are given 

in Table |Vll 

From Fig. [2] we can see that adding a pulsar compo- 
nent can roughly reproduce the AMS-02 positron frac- 
tion data and the Fermi/HESS total electron spectra. 
However, it also shows a tension between the AMS-02 
positron fraction and the Fermi/HESS electron spectra, 
under the current framework. Quantitatively, the best 
fitting x^ is about 279 for fit I-a. For 151 dof such a x^ 



means ^ 6.1a deviation from what expected. The best 
fitting model is a little over-estimated compared with the 
AMS-02 data, while under-estimated compared with the 
Fermi data. This is clearly seen from the fits Il-a, in 
which the constraints of Fermi/HESS data arc removed. 
We sec that the fit Il-a gives good fit to the AMS-02 data, 
but deviate from the Fermi data obviously. 

Compared with the fitting results with PAMELA 
positron fraction data |j|, the spectrum of the pulsar 
component becomes much softer (with power-law index 
^ 1.9), which may be more reasonable according to the 
pulsar modeling [7|. The contribution of e from the 
pulsars is also smaller than previous estimated according 
to the PAMELA data. Our fit shows that up to TeV the 
positron fraction is only ~ 20%, while it is more than 
30% or even reaching 40% according to the fitting to the 
PAMELA data. 

Figs. [3] and m give the results for the DM annihilation 
scenario. The tension between AMS-02 and Fcrmi/HESS 
data becomes stronger and the fits get worse than the 
pulsar scenario. For DM annihilation into fi^fi", the 
reduced x^ for fit I-a is as high as 3.3. The reason of 
the poor fit is that the positron spectrum from DM an- 
nihilation to a pair of /i^/i^ is too hard. This can be 
seen from the top-left panel of Fig. [3l The DM com- 
ponent will over-produce high energy positrons (> 100 
GeV) but under-producc positrons at tens of GcV. For 
fits I-a, heavy DM with mass '--'2 — 3 TeV is required due 
to the constraints of Fermi/HESS data. If we throw away 
the Fermi/HESS e± data (fit Il-a), we find the AMS-02 
data tend to favor lighter DM particles. In this case the 
10 - 50 GeV AMS-02 data can be better reproduced, 
however, the data above 100 GeV are still difficult to be 
explained. 

The situation for t'^t~ final state is better. Simi- 
lar with the ^i^ n~ channel, the positron spectrum from 
tauon decay is still too hard. Relaxing the constraints 
from Fcrmi/HESS data we find a lighter DM mass, 
- 550 - 750 GeV is favored. In this case the AMS-02 
data can be fitted relatively well. It is interesting to note 
that even only the AMS-02 data and PAMELA electron 
data are considered (fit Il-a), the mass of DM particles 
can be constrained in a small region. Such a strong con- 
straint comes from the very high precise AMS-02 data at 
lower energy (a few tens GeV). Therefore, we should not 
take these values too seriously. The change of the back- 
ground, such as taking different propagation parameters 
or hadronic interaction models (we will discuss this issue 
later), may affect the detailed parameters. 

Fig. [5] gives the la and 2a contour for the DM mass 
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FIG. 2: The positron fraction (upper) and electron spectra (lower) for the background together with a pulsar-like component 
of the exotic e . The panels from left to right are for fits I-a and Il-a respectively. References of the data: positron fraction 
— AMS [5l^EAT94+95 [H, HEATOO [13, PAMELA 0, AMS-02[i| ; electron — PAMELA fll, ATIC ^, HESS ^^, 
Formi-LAT mil. 



TABLE HL Fitting results of pulsar-like model with proton spectrum fixed 





I-a 


Il-a 




best 


mean 


best 


mean 


l0g(Ae'") 


-8.978 


-8.974 ± 0.005 


-8.925 


-8.921 ±0.011 


71 


1.504 


1.512 ±0.010 


1.708 


1.704 ±0.084 


72 


2.645 


2.652 ±0.010 


2.794 


2.796 ± 0.028 


log(pC7MeV) 


3.599 


3.587 ± 0.022 


3.597 


3.600 ± 0.046 


log{A^J) 


-24.867 


-24.918 ±0.146 


-25.257 


-25.226 ±0.562 
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1.912 


1.903 ±0.029 


1.856 


1.863 ±0.116 


log(p,/MeV) 


6.640 


6.632 ±0.111 


5.927 


6.097 ± 0.412 


Ce + 


1.272 


1.327 ± 0.075 


2.206 


2.222 ± 0.242 


<t)/MY 


500 


527 ± 30 


818 


830 ± 72 



"Normalization at 25 GeV in unit of cm s sr" 



Normalization at 1 MeV in unit of cm s ^MeV 



^MeV^ 
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and annihilation cross section. But we should keep in 
mind that such results should not be considered statisti- 
cally meaningful as the fits are quite bad. The solid lines 
shown in Fig. [S] are the exclusion limits derived by the 
Fermi 7-ray observations of the Galactic center [53| and 
dwarf galaxies [60| • We can see that 7-rays tend to give 
strong constraints on the DM scenario, especially for the 
T'^T" final state. 



We further note that for the DM scenario, the param- 
eter <f) is very large. The solar modulation potential is 
assumed to vary between 300 and 1000 MV in these fits. 
From Tables IIIII - |V] we see that almost in all cases the 
modulation potential tends to the upper end. This might 
be inconsistent with the fact that PAMELA and AMS-02 
work approaching the solar minimum. 

Since there might be discrepancy between the AMS- 
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FIG. 3: Same as Fig. [2] but the exotic e are assumed to be from DM annihilation. The annihilation channel is fJ,'^ fJ. 
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FIG. 4: Same as Fig. [2] but the exotic e are assumed to be from DM annihilation. The annihilation channel is t^t 



02 and Fermi/HESS data, we give the expected fluxes of pure positrons and electrons seperately in Fig. [S]based on 



TABLE IV: Fitting results of DM annihilation to /i^^ with proton spectrum fixed 





I-a 


Il-a 




best 


mean 


best 


mean 


log(A.") 


-8.916 


-8.916 ±0.003 


-8.915 


-8.918 ±0.006 


71 


1.894 


1.870 ±0.036 


1.896 


1.868 ±0.036 


72 


2.839 


2.839 ± 0.006 


2.904 


2.906 ±0.014 


log(pL/MeV) 


3.608 


3.592 ± 0.035 


3.692 


3.674 ± 0.037 


log(m;,/GeV) 


3.371 


3.368 ± 0.039 


2.415 


2.423 ± 0.039 


log((a-w)/cm'*s^"'") 


-22.307 


-22.313 ±0.067 


-24.169 


-24.166 ± 0.069 


Ce + 


2.881 


2.881 ± 0.030 


3.052 


3.038 ± 0.047 


4,/MV 


999 


996 ±4 


999 


991 ±8 



"Normalization at 25 GeV in unit of cm ^s ^sr ^MeV ^. 



TABLE V: Fitting results of DM annihilation to t^t with proton spectrum fixed 
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mean 


best 


mean 


log(A.") 


-8.915 


-8.916 ±0.003 


-8.907 


-8.909 ± 0.006 


71 


1.879 


1.878 ± 0.037 


1.869 


1.817 ±0.070 


72 


2.813 


2.813 ±0.007 


2.863 


2.856 ±0.015 


log(pL/MeV) 


3.570 


3.571 ±0.037 


3.637 


3.608 ± 0.058 


log(m;,/GeV) 


3.667 


3.665 ± 0.045 
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log({cr^)/cm^s"^) 


-21.699 


-21.703 ± 0.073 


-23.261 


-23.269 ± 0.064 


Ce + 


2.773 


2.769 ± 0.035 
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2.844 ± 0.068 


(t)/M\f 


999 


994 ±5 


998 


974 ± 20 



"Normalization at 25 GeV in unit of cm 
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FIG. 5: Icr and 2cr confidence regions on the DM mass and cross section plane, for the fits I-a and Il-a respectively. The left 
panel is for fJ,'^ i-i~ channel, and the right panel is for t'^t~ channel. The solid lines show the 95% upper limit of Fermi 7-ray 
observations of the Galactic center (with normalization of the local density corrected) [591 ] and dwarf galaxies [601 ]. 



the fit Il-a (i.e., without Fermi/HESS data). It is shown 
that the electron spectrum shows an almost featureless 
behavior and the extra component is difficult to be iden- 
tified from the pure electron spectrum. It is more promi- 
nent to see the extra component in the positron spec- 



trum, which has a very hard spectrum (~ E~'^'^) above 
^10 GeV. The following measurement of the positron 
and electron fluxes by AMS-02 can directly test such re- 
sults. 



TABLE VI: Summary of fitting x^/dof. Note that the AMS- 
02 data above 5 GeV and PAMELA proton data below 150 
GeV are used to calculate the x^- 





pulsar 


DM ifi+^i-) 


DM (r+T-) 


La 
ILa 


278.7/151 
51.5/80 


506.7/152 
83.1/81 


496.5/152 
56.7/81 


Lb 

ILb 


288.0/205 
83.0/134 


615.3/206 
238.7/135 


584.6/206 
164.3/135 




FIG. 6: Expected fluxes of electrons (solid) and positrons 
(dashed) for the fit ILa. The higher one in each group is 
the local interstellar (LIS) flux and the lower one represents 
the flux after the solar modulation. Grey points show the 
PAMELA measurement of electrons [ll| and CAPRICE mea- 
surement of positrons [63| . 



positron spectrum is too hard. The la and 2a confidence 
level contours on the m^ — (av) plane for the fits are 
shown in Fig. 1101 The parameter regions differ only 
slightly from that derived in the previous subsection (Fig. 
[5]) . The strong constraints on the DM model from 7-rays 
are not changed. 

From Figs. |8]and|9]we note that the proton spectrum 
can not be well fitted for the DM scenario. This is also 
due to the hard positron spectrum from DM annihila- 
tion into muons and tauons. A harder proton spectrum 
will produce more positrons above ~ 10 GeV, which will 
compensate the lack of positrons from the hard spectrum 
of the DM component. If we reduce the constraints from 
Fermi/HESS data (fits Il-b), we see that the proton spec- 
trum fits the data better. 

We also note from Figs. |8] and [9] that when not includ- 
ing Fermi/HESS data the DM scenario docs not give a 
good description to the high energy end of the AMS-02 
data, which is different from the pulsar case shown in Fig. 
[71 This is because the positron spectrum is determined 
once the mass of DM and its annihilation final states are 
given. To have a minimum x^ the mass of DM (i.e. the 
shape of the positron spectrum) is usually determined by 
the data at tens of GeV where the errors are very small, 
instead of the behavior of the high energy end data. 

The solar modulation potential and Ce+ at the DM 
case are larger than that in the pulsar case. The reason is 
that the DM spectrum is always harder than the pulsar 
case. To fit the AMS-02 data at tens of GeV, which are 
very precise, larger (j) and Ce+ can give relatively higher 
fiux of positrons in this energy range. 



B. Relaxing the proton spectrum 



Further tests 



Since the proton spectrum will affect the secondary 
positron production, and also the determination of the 
solar modulation parameter, we take the proton spec- 
trum into account and redo the fits (labeled as I-b and 
Il-b). The fitting results with the best fitting parameters 
are presented in Figs. [7]l9l The best fitting results, mean 
values and the la uncertainties of the model parameters 
are listed in Tables IVIHIXl and the best fitting x^ over 
dof are also presented in Table IVII Qualitatively we find 
that the results are similar with that in the previous sub- 
section. 

For the pulsar scenario, the minimum x^ ioi fit I-b 
is about 288, which corresponds to a ^ 3.9a deviation 
from what expected for 205 dof. Compared with the 
fits of fixed proton spectrum, the injection spectrum of 
positrons from pulsars is a little bit softer here. This 
is because the proton spectrum here is also softer than 
that in the previous subsection. Therefore a softer pulsar- 
induced positron spectrum is required to give more tens 
of GeV positrons and to compensate the effect of a softer 
proton spectrum. 

The DM scenario fits worse than the pulsar scenario. 
As we have discussed, the reason is that the DM-induced 



To better understand how soft a positron spectrum 
from the extra sources is needed, we show in Fig. [11] the 
2(7 range of source spectra from pulsars at the solar loca- 
tion with shaded regions. The 2a range is defined with 



Ax' 



X 



/Vniin 



22.7 for 13 fitting parameters. For 



comparison the DM induced positron spectra for fJ,'^ ^~ 
T~^T~ , W'^W' and hh channels arc also shown. The mass 
of DM particle is taken to be m^ = 1 TeV and a free flux 
normalization is adopted. It is shown that the positron 
spectra from DM annihilation in the /i^/i^ and t^t~ 
channels are much harder than the pulsar component. 
For the W^W~ and hh channels the spectrum is softer 
and we would expect a better fit to the data. 

As a test we run the fit Il-a with DM annihilation 
to W~^W~ and hh final states. We find the x^ values 
become slightly smaller (^ 52.7 for both W~^W~ and bb) 
than that of t'^t~. The positron fraction and electron 
spectrum for the best fitting parameters are shown in 
Fig. [T21 It is shown that the AMS-02 positron fraction 
data can be reproduced in this case. 

We know that the PAMELA antiproton data and 
Fermi 7-ray data set very stringent constraints on the 
DM annihilation into quark and gauge boson final 




galprop 

AMS 

BESS 

ATIC2 

PAMELA 

CREAM 




E^ (GeV/nucleon) 



E^ (GeV/nucleon) 




FIG. 7: From top to bottom: the proton, positron fraction and electron spectra for the background together with a pulsar-hke 
component of the exotic e*. The left and right panels are for fits I-b and Il-b respectively. 



states[43, H, Ifi^ilfiJ]- Fig. [T3l shows the two dimensional 
contours on the m^ — (av) plane for the W~^W~ and bb 
channels and the 95% exclusion limits on DM annihila- 
tion to bb and W~^W~ channels by Fermi observation of 
dwarf galaxies [601 ■ We sec that the current positron frac- 
tion data of AMS-02, together with other e^, antiproton 
and 7-ray data, give strong constraints on the present 
DM scenario to explain the e^ excesses. 



V. DISCUSSION 
A. Uncertainties of the theoretical model 

Since the AMS-02 data are very precise, any uncertain- 
tics previously thought to be not important may affect 
the fitting result. As we note above the positron spec- 
trum around tens of GeV may affect the fitting result 
sensitively. Therefore the change of the shape of back- 
ground may affect the fitting results. There are quite a 
few sources of such uncertainties, such as the propaga- 
tion parameters, the hadronic interaction models, and so 
on. 

In this work we use the Kamae et al. (2006) parame- 
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TABLE VII: Fitting results of pulsar-like model with proton spectrum relaxed 





I-b 


Il-b 




best 


mean 


best 


mean 


log(A/) 


-8.323 


-8.327 ± 0.005 


-8.328 


-8.330 ± 0.006 


yi 


1.789 


1.797 ±0.019 


1.882 


1.885 ±0.023 


Vl 


2.378 


2.388 ±0.011 


2.410 


2.415 ±0.019 


logKr/MeV) 


4.040 


4.064 ± 0.029 


4.099 


4.121 ±0.031 


log(A,') 


-8.977 


-8.979 ± 0.004 


-8.941 


-8.948 ± 0.008 


71 


1.504 


1.505 ± 0.004 


1.535 


1.554 ± 0.040 


72 


2.647 


2.645 ±0.011 


2.720 


2.722 ±0.018 


log(pL/MeV) 


3.615 


3.614 ±0.015 


3.614 


3.631 ±0.016 


log(Ap,,^) 


-25.104 


-25.012 ±0.132 


-24.494 


-24.411 ±0.335 


a 


1.864 


1.881 ±0.026 


1.985 


2.006 ± 0.068 


log(p,/MeV) 


6.512 


6.562 ± 0.093 


6.205 


6.442 ± 0.285 


Ce + 


1.306 


1.276 ±0.062 


1.503 


1.468 ±0.127 


(^/MV 


490 


489 ± 21 


623 


614 ±37 



"Normalization at 100 GeV in unit of cm s 



Normalization at 25 GeV in unit of cm s ^sr 



sr^^MeV"^ 
^MeV"\ 



''Normalization at 1 MeV in unit of cm ^s ^MeV 



TABLE VIII: Fitting results of DM annihilation into yt^ ^i with proton spectrum relaxed 





Lb 


ILb 




best 


mean 


best 


mean 


log(A/) 


-8.275 


-8.276 ± 0.004 


-8.299 


-8.298 ± 0.004 


V\ 


1.936 


1.947 ±0.031 


2.037 


2.026 ± 0.022 


Vl 


2.271 


2.281 ± 0.009 


2.364 


2.365 ±0.010 


log(p^,/MeV) 


3.919 


3.929 ± 0.045 


4.102 


4.101 ±0.035 


log(A,') 


-8.953 


-8.950 ± 0.006 


-8.938 


-8.944 ± 0.008 


71 


1.555 


1.588 ± 0.054 


1.789 


1.754 ± 0.046 


72 


2.768 


2.773 ±0.011 


2.890 


2.895 ±0.016 


log(p^,/MeV) 


3.578 


3.573 ± 0.027 


3.679 


3.670 ± 0.026 


log(m;,/GeV) 


3.330 


3.338 ± 0.045 


2.390 


2.417 ±0.038 


log((crii)/cm'*s^"'") 


-22.397 


-22.381 ±0.076 


-24.142 


-24.117 ±0.060 


Ce+ 


1.996 


2.037 ± 0.064 


2.541 


2.527 ±0.070 


,^/MV 


702 


729 ± 39 


879 


866 ± 36 



"Normalization at 100 GeV in unit of cm"^s"^sr"^MeV"^ 
''Normalization at 25 GeV in unit of cm~'^s~^sr~^MeV~^. 



terization of the -p-p collision [37[ . As shown in [6J| there 
were remarkable differences between different hadronic 
models. The Kamae et al. (2006) parameterization in- 
cluded more processes than before and was calibrated 
with recent data [37[. However, it depends strongly on 
the Monte Carlo simulations. Therefore we also test the 
fitting with old "pp collision parameterization |65| . The re- 
sults show quantitative difference from that shown above. 
But, all the conclusions made in the previous section keep 
unchanged. 



Additional spectral break of the primary 
electrons 



As revealed by ATIC, CREAM and PAMELA mea- 
surements, the nuclei spectra have a hardening above 
- 200 GV [H m [13. It is possible that the primary 
electron spectrum also has a similar hardening at high 
energies j66| . By including such a modification of the pri- 
mary electron spectrum may soften the tension between 
the AMS-02 and Fermi data. This will be investigated in 
future works. 
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FIG. 8: Same as Fig. [7] but the exotic e are assumed to be from DM annihilation. The annihilation channel is fJ,'^ fJ. 



C. Two or more components of the extra e 
sources 

It is possible that there are more than one compo- 
nents of the primary sources of positrons. For example 
for the pulsar scenario, the far away pulsar population 
may contribute to a "background" component, and sev- 
eral nearby pulsars may give very distinct contributions 
to the positron spectrum [5|. Such a picture may also 
help to eliminate the tension between AMS-02 data and 
Fermi/HESS e^ data. The "background" pulsars may 
give most contributions to the AMS-02 positron excesses, 
while the nearby sources can contribute mainly at high 
energies to reproduce the Fermi/HESS data. We leave 
the detailed modeling in future works. 



D. To improve the DIVI scenario 

Finally we discuss the possibilities to improve the DM 
scenario to give better explanation of the current data. 
First we need softer and broader spectrum of positrons 
from DM. Therefore if the annihilation final state is not 
two-body state but four-body, eight-body etc., softer 
positron spectrum may be generated. A mixture of Icp- 
tonic channels and hadronic channels may also give a 
broader spectrum (see Fig. Ilip . Second we have to 
consider how to avoid the constraints from 7-rays and 
antiprotons. Decaying DM scenario is better since the 
constraints from 7-rays are less stringent [63, [G^. As 
for the antiproton constraints, very massive DM particle 
may avoid the current bounds set by the PAMELA data 
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FIG. 9: Same as Fig. [7] but the exotic e are assumed to be from DM annihilation. The annihilation channel is t^t 



m- 



VI. SUMMARY 



In summary in this work we give a systematical inves- 
tigation of the models to explain the cosmic e excesses, 
based on the newest AMS-02 data of the positron fraction 
and other data from PAMELA, Fermi and HESS. Both 
the pulsar scenario and DM scenario as the extra primary 
e sources are studied. Our findings are as follows. 



positrons from the extra sources than that needed 
by Fermi/HESS data. 

• Pulsar models can roughly fit the data. The spec- 
tral index of the positrons injected by pulsars is 
about E^^, which is much softer than that derived 
when using the PAMELA positron fraction data. 

• If we fit only the AMS-02 positron fraction and 
PAMELA electron spectrum data both the pul- 
sar and DM annihilates/decays into t^t~ , W~^W~ 
and bb can fit the data. 



• In the present framework there is tension be- 
tween the AMS-02 data and Fermi/HESS data. 
The AMS-02 positron fraction data requires less 



• Due to the strong constraints from antiprotons 
and 7-rays, the DM annihilation scenario (into two 
body final states) seems not easy to be consistent 
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TABLE IX: Fitting results of DM annihilation into t^t with proton spectrum relaxed 





I-b 


Il-b 




best 


mean 


best 


mean 


log(A/) 


-8.278 


-8.282 ± 0.004 


-8.310 


-8.310 ±0.005 


Vl 


1.900 


1.892 ±0.023 


1.962 


1.951 ±0.030 


Vl 


2.281 


2.285 ± 0.008 


2.387 


2.381 ±0.012 


log(pL/MeV) 


3.936 


3.920 ± 0.032 


4.136 


4.101 ±0.037 


log(A,^) 


-8.962 


-8.961 ±0.004 


-8.940 


-8.944 ± 0.008 


71 


1.512 


1.534 ±0.026 


1.637 


1.604 ± 0.056 


72 


2.730 


2.731 ±0.011 


2.794 


2.792 ± 0.021 


log(p§,/MeV) 


3.558 


3.569 ± 0.020 


3.633 


3.610 ±0.027 


log(m;,/GeV) 


3.555 


3.522 ± 0.075 


2.657 


2.670 ± 0.031 


log((cr7;)/cni^s"^) 


-21.884 


-21.929 ±0.116 


-23.222 


-23.216 ±0.048 


Ce + 


1.856 


1.875 ±0.051 


2.167 


2.151 ±0.094 


(^/MV 


650 


649 ± 29 


736 


733 ± 47 



"Normalization at 100 GeV in unit of cm"^s"^sr"^MeV" 
'Normalization at 25 GeV in unit of cm~^s~^sr~^MeV~ 
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FIG. 10: lo" and la confidence regions on the DM mass and cross section plane, for the fits I-b and Il-b respectively. The left 
panel is for fj.'^ jJ.~ channel, and the right panel is for t'^t~ channel. The solid lines show the 95% upper limit of Fermi 7-ray 
observations of the Galactic center (with normalization of the local density corrected) [591 ] and dwarf galaxies pU^ . 



with all of the current data. 

Our study illustrates the remarkable potential on un- 
derstanding the physics of the e^ excesses from the very 
precise measurement done by AMS-02. According to the 
PAMELA positron fraction data, the pulsar and DM sce- 
narios are almost identical Q . Given the AMS-02 precise 
measurement, the differences between different scenarios 
become to appear. We are looking forward to more data 
with higher precision from AMS-02 to further shed light 
on the understanding of the fundamental questions of 



both astrophysics and particle physics. 
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